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[l]  Satellite  observations  of  polar  stratospheric  clouds 
(PSCs)  over  Antarctica  in  June  2003  revealed  small  nitric 
acid  trihydrate  (NAT)  particles  forming  suddenly  along  the 
vortex  edge.  Models  suggest  the  trigger  was  mountain  waves 
over  the  Antarctic  Peninsula  (AP)  forming  ice  for  NAT 
nucleation.  We  test  this  hypothesis  by  analyzing  perturbations 
in  stratospheric  radiances  from  the  Atmospheric  Infrared 
Sounder  (AIRS).  AIRS  data  show  mountain  waves  over  the 
AP  on  10-14  June,  with  no  resolved  wave  activity  before  or 
after.  Peak  wave  temperature  amplitudes  derived  from 
independent  40  hPa  channels  all  return  values  of  10-12  K, 
in  agreement  with  values  used  to  model  this  NAT  event. 
These  observations  support  a  NAT  wake  from  a  small  region 
of  mountain  wave  activity  over  the  AP  as  the  source  of 
this  circumpolar  NAT  outbreak.  Citation:  Eckermann,  S.  D., 
L.  Hoffmann,  M.  Hopfner,  D.  L.  Wu,  and  M.  J.  Alexander  (2009), 
Antarctic  NAT  PSC  belt  of  June  2003:  Observational  validation  of 
the  mountain  wave  seeding  hypothesis,  Geophys.  Res.  Lett.,  36, 
L02807,  doi:  1 0. 1 029/2008GL036629. 

1.  Introduction 

[2]  Each  winter,  PSCs  form  as  polar  stratospheric  temper¬ 
atures  drop  below  ~195  K.  Three  nominal  PSC  types  - 
nitric  acid  trihydrate  (NAT),  supercooled  ternary  solutions 
(STS)  and  ice  -  can  exist  at  progressively  lower  temper¬ 
atures  T  ( Tice  <  TSts  <  Tnat).  The  cold  interior  of  the 
Antarctic  polar  vortex  supports  a  persistent  polar  cap  of 
PSCs  that  converts  reservoir  chlorine  species  into  activated 
forms  that  destroy  ozone.  A  stable  meridional  temperature 
gradient  drives  strong  vortex-edge  winds  that  limit  lateral 
mixing. 

[3]  This  temperature  gradient  can  cause  PSC  types  to 
vary  with  equivalent  latitude.  Near  the  vortex  edge  where 
Tsts  <  T  <  Tnat,  a  circumpolar  “freezing  belt”  of  NAT 
PSCs  can  exist  [Tabazadeh  et  al.,  2001].  If  NAT  particles 
form  then  remain  within  freezing-belt  latitudes  for  several 
or  more  days,  they  can  grow  to  sizes  large  enough  to 
sediment  to  lower  sublimation  altitudes.  Over  time  this 
process  can  denitrify  the  stratosphere,  accentuating  spring¬ 
time  ozone  loss  by  removing  odd  nitrogen  compounds  that 
sequester  activated  chlorine.  How,  when  or  if  NAT  forms  in 
these  regions  depends  on  dominant  nucleation  mechanisms, 
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which  remain  uncertain  and  controversial  [Lowe  and 
MacKenzie,  2008]. 

2.  June  2003  NAT  PSC  Outbreak 

[4]  The  Michelson  Interferometer  for  Passive  Atmo¬ 
spheric  Sounding  (MIPAS)  performed  high  spectral  resolu¬ 
tion  limb  measurements  of  the  Antarctic  winter  stratosphere 
in  2003.  A  color-ratio  algorithm  described  by  Hopfner  et  al. 
[2006a]  enabled  PSCs  of  different  types  to  be  discriminated. 

[5]  Figure  1,  adapted  from  Hopfner  et  al.  [2006b],  maps 
all  PSCs  measured  by  MIPAS  on  6  selected  days  in  June 
2003.  Red  squares  show  solid  NAT  PSCs  of  <3  pm  radius. 
A  NAT  freezing  belt  could  exist  throughout  early  June 
(orange  band),  yet  no  NAT  PSCs  were  observed.  This 
suddenly  changed  on  10-11  June  when  small  NAT  particles 
streamed  from  a  genesis  region  near  the  Antarctic  Peninsula 
(AP)  (Figure  Id),  a  feature  observed  more  strongly  on  12  June 
and  leading  3-4  days  later  to  a  freezing  belt  of  solid  NAT 
particles  most  of  the  way  around  Antarctica  (Figure  If). 

[6]  Hopfner  et  al.  [2006b]  modeled  this  event  using  a 
stratospheric  trajectory  model  with  microphysics  under 
various  proposed  NAT  nucleation  mechanisms.  No  simula¬ 
tion  reproduced  these  observations  until  mountain  wave- 
induced  temperature  perturbations  were  considered.  A 
mountain  wave  forecasting  algorithm  [Eckermann  et  al., 
2006a]  predicted  stratospheric  temperature  perturbations  of 
~±10  K  over  the  AP  on  10-14  June  with  little  or  no  wave 
activity  before  and  after.  By  inserting  these  perturbations 
into  their  model,  the  observations  were  well  reproduced  via 
NAT  nucleating  heterogeneously  on  ice  formed  in  the 
cooling  phases  of  AP  mountain  waves.  The  simulation 
reproduced  the  sudden  onset  of  NAT  on  11  June  near  the 
AP.  The  persistence  of  wave  activity  through  14-15  June 
provided  a  constant  source  of  ice  on  which  NAT  could 
nucleate,  and  a  continuous  wake  of  solid  NAT  that  swept 
out  a  circumpolar  NAT  freezing  belt  by  15  June. 

[7]  This  explanation  relies  on  model  predictions  of  10- 
15  K  stratospheric  mountain  wave  amplitudes  over  the  AP 
during  this  event.  Observational  validation  of  these  moun¬ 
tain  wave  predictions  remains  the  key  missing  test  of  this 
hypothesis. 

3.  Mountain  Waves  in  AIRS  Radiances 

[8]  AIRS  scans  cross-track  in  90  step-and-stare  measure¬ 
ments  distributed  symmetrically  either  side  of  nadir.  Hori¬ 
zontal  measurement  footprints  are  ^13  x  13  km2  at  nadir 
and  ~3 1  x  21  km2  at  the  far  off-nadir  views,  with  a  total 
cross-track  swath  of  ^1650  km.  As  scans  repeat,  “pushb- 
room”  radiance  images  are  swept  out.  Despite  having  deep 
vertical  weighting  functions,  the  high  horizontal  resolution 
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Figure  1.  MIPAS  Antarctic  PSC  measurements  on 
indicated  dates  (adapted  from  Hopfner  et  al.  [2006b])  at 
~21  km  altitude.  Red  squares  show  NAT  particles  with 
radii  <3  pm,  blue  triangles  show  ice,  green  circles  show 
other  PSC  types  (either  STS,  large  NAT  or  thin  ice).  Red, 
green  and  blue  curves  mark  NAT,  STS  and  ice  equilibrium 
formation  thresholds,  respectively,  based  on  analysis 
temperatures.  Orange  band  marks  the  region  Tsrs  <  T  < 
Tnat- 

and  radiometric  precision  of  these  measurements  allow 
short  deep  gravity  waves  to  be  resolved,  particularly  in 
stratospheric  channels  negligibly  impacted  by  tropospheric 
moisture  [Alexander  and  Barnet,  2007;  Alexander  and 
Teitelbaum,  2007;  Eckermann  et  al.,  2007]. 


[9]  Since  AIRS  has  over  2300  channels,  to  both  thin  the 
data  and  reduce  noise  we  averaged  Level  lb  radiances 
from  a  subset  of  thermal  channels  with  similar  vertical 
weighting  functions.  We  focus  initially  on  an  80  hPa 
product  derived  by  averaging  radiances  from  channels 
92,  98,  104,  105,  110,  111,  116,  117,  122,  123,  128, 
129,  134,  and  140.  Radiances  with  horizontal  scales 
>600  km  were  fitted  as  in  Eckermann  et  al.  [2006b]  then 
subtracted  from  the  data  to  leave  small-scale  residuals.  We 
applied  3x3  point  smoothing  to  reduce  noise  further  and 
accentuate  any  geophysical  wave  signals. 

[10]  Figure  2  plots  a  time  sequence  of  80  hPa  radiance 
residuals  overthe  AP  for  overpasses  from  8- 15  June  2003.  On 
8-9  June  2003  there  is  no  resolved  wave  structure.  Late  on 
10  June  we  see  first  evidence  of  weak  wave  oscillations  in 
Figure  2d.  Oscillations  intensify  and  spread  on  11  June 
(panels  e-f),  peak  in  amplitude  on  12-13  June  (panels  i-k), 
weaken  on  14  June  and  disappear  by  15  June.  This  time 
evolution  is  qualitatively  consistent  with  the  mountain  wave 
forecasts  in  Figure  3  of  Hopfner  et  al.  [2006b]. 

[11]  Similar  maps  of  brightness  temperature  perturbations 
T B  from  channel  72,  which  peaks  ~6  km  higher  near 
30  hPa,  revealed  clear  oscillations  (not  shown)  that  are 
~180°  out  of  phase  with  those  in  Figure  2,  indicating  a 
vertical  wavelength  A-  «  12  km.  A  stationary  hydrostatic 
mountain  wave  has  a  theoretical  vertical  wavelength  Xz  = 
2ttU/N,  where  U  is  horizontal  wind  speed  orthogonal  to 
wave  phase  fronts  and  N  is  buoyancy  frequency.  On  12- 
13  June,  Met  Office  analyses  over  the  AP  at  30-80  hPa 
yield  U  «  35  m  s-1  and  N  «  0.018  rad  s-1,  and  thus  Az  « 
12  km,  consistent  with  the  observational  estimate.  Thus  the 


[a) 8  June  2003 0418 UTC:  80.0 hPa  [b]  9  June 2003  1824 UTC: 80.0  hPa  [c]  10 June  2003  0406  UTC:  80.0 hPa (d)  10 June 2003  1907  UTC:  80.0  hPa[e]  11  June 20030310 UTC:  80.0 hPa 
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Figure  2.  Brightness  temperature  perturbations  TB,  extracted  from  80  hPa  AIRS  radiances  (see  text  for  details),  for  select 
overpasses  of  the  AP,  in  chronological  order  from  (a)  top-left,  8  June  0418  UTC,  to  (o)  bottom-right,  15  June  0424  UTC. 
Date  and  time  of  each  overpass  is  given  in  plot  title,  as  well  as  maximum  and  minimum  fB  values. 
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Figure  3.  (a)  Gray  curves  show  vertical  weighting 
functions  for  AIRS  channels  64,  88,  90,  94,  100,  106  and 
118,  spanning  wavenumbers  665.015-678.839  cm-1,  and 
the  corresponding  mean  (black  broken  curve).  Solid  curve 
shows  weighting  function  for  channel  71  (666.773  cm-1), 
(b)  Corresponding  attenuations  e  as  a  function  of  Az. 

observed  oscillations  are  consistent  with  stationary  AP 
mountain  waves. 

4.  Wave  Temperature  Amplitudes  at  PSC  Height 

[12]  To  assess  the  impact  of  these  mountain  waves  on 
NAT  formation,  we  must  infer  peak  wave  temperature 
amplitudes  T  at  the  ~2 1  km  geometric  height  (~40  hPa) 
of  the  NAT  PSCs  in  Figure  1. 

[13]  Gray  curves  in  Figure  3a  plot  vertical  nadir  weight¬ 
ing  functions  of  a  collection  of  AIRS  1 5  pm  channels  that 
peak  near  40  hPa  with  similar  properties.  We  average 
radiances  from  these  channels,  as  for  the  80  hPa  radiances 
in  Figure  2,  into  a  mean  radiance  scene  with  reduced  noise. 
As  a  cross  check,  we  separately  analyze  radiances  from 
channel  71  that  peaks  more  sharply  at  40  hPa.  These 
weighting  functions  were  computed  using  the  radiative 
transfer  model  of  L.  Hoffmann  and  M.  J.  Alexander 
(Retrieval  of  stratospheric  temperatures  from  AIRS  radiance 
measurements  for  gravity  wave  studies,  submitted  to  Jour¬ 
nal  of  Geophysical  Research,  2008)  using  a  high-latitude 
winter  temperature  profile. 

[14]  These  broad  weighting  functions  attenuate  the  ampli¬ 
tude  of  any  gravity  wave  in  the  scene  as 

Yb  =  <Xz)T  (1) 

where  T B  is  (measured)  peak  brightness  temperature 
amplitude,  e  is  the  attenuation  and  overbars  denote  back¬ 
ground  values.  Figure  3b  plots  e(Az)  for  these  40  hPa 
channels  by  scaling  the  relative  radiance  response  functions 
of  Hoffmann  and  Alexander  (submitted  manuscript,  2008) 
using  equation  (5)  of  Alexander  and  Barnet  [2007]  and  T  k, 
Tbk.  195  K.  For  X-  ~  12  km,  t B  is  ~7%  of  the  actual  wave 
temperature  amplitude  T  for  the  multichannel  40  hPa 
radiances,  and  ~17%  for  channel  71  radiances. 

[15]  Figure  4  plots  40  hPa  7y  perturbations,  derived  by 
scaling  TB  fields  by  e-1  values  at  Az  =  12  km,  for  over¬ 
passes  on  12  and  13  June.  Left  panels  used  multichannel 
40  hPa  radiances,  right  panels  used  channel  71  radiances. 
Despite  very  different  T B  and  e  values,  each  returns  the 
same  independent  estimates  of  T  x;  1 2  K  on  1 2  June  and 
«10  K  on  13  June.  For  an  error  of  ±1  km  in  Az,  the 


corresponding  uncertainty  in  e  (see  Figure  3b)  leads  to  error 
bars  of  ~±2-4  K  in  T . 

[16]  These  mean  40  hPa  T  estimates  are  probably  lower 
bounds.  The  finite  horizontal  width  of  individual  AIRS 
measurements  reduces  wave  amplitude.  More  importantly, 
our  3x3  point  smoothing,  while  reducing  small-scale 
noise,  also  reduces  TB  values.  Comparisons  of  smoothed 
and  unsmoothed  scenes  (not  shown)  suggest  reductions  of 
~20-30%.  Thus  actual  values  could  be  nearer  15  K. 

5.  Discussion  and  Summary 

[17]  Our  inferred  40  hPa  amplitudes  of  10-15  K  for 
the  AP  mountain  waves  measured  by  AIRS  validate  the 
t  =  10-15  K  values  used  by  Hopfner  et  al.  [2006b], 
based  on  hindcasts  with  the  Mountain  Wave  Forecast 
Model  [Eckermann  et  al.,  2006a].  When  inserted  into 
their  microphysical  model,  Hopfner  et  al.  [2006b]  showed 
that  AP  waves  of  this  amplitude  formed  ice.  Solid  NAT 
particles  nucleated  on  this  ice  and  survived  downstream. 
Only  this  nucleation  scenario  reproduced  the  observations 
in  Figure  1.  Thus  these  AIRS  observations  validate  the 
Hopfner  et  al.  [2006b]  mountain  wave  explanation  of  this 
observed  NAT  outbreak. 

[is]  This  event  represents  a  remarkable  observational 
verification  of  the  viability  of  the  “mother  cloud” 
[Fueglistaler  et  al.,  2002]  or  “salt  shaker”  [Dhaniyala  et 
al.,  2002]  theory  of  vortex-wide  solid  NAT  formation  and 
growth  driven  by  ice  formed  in  isolated  pockets  of  strato¬ 
spheric  mountain  wave  activity.  The  AP  mountain  waves 
observed  here  have  horizontal  wavelengths  A/,  «  300  km. 
The  cold  phases  where  ice  forms  are  perhaps  no  more  than 
30  x  100  km2  in  horizontal  extent.  Yet,  in  this  case,  these 
small  regions  of  mountain  wave-induced  ice  production 
prove  sufficient  to  convert  the  edge  of  the  Antarctic  vortex 
from  totally  NAT  free  (Figures  la-lc)  to  being  fully 
populated  with  NAT  PSCs  (Figures  le—  If)  over  a  period 
of  less  than  a  week. 

[19]  Almost  all  research  on  mountain  wave-induced  NAT 
has  focused  on  the  Arctic,  where  temperatures  are  wanner 


Figure  4.  Inferred  40  hPa  T  perturbations  at  1855  UTC  on 
12  June  (top  row)  and  1938  UTC  on  13  June  (bottom  row), 
derived  by  scaling  T B  fields  by  e-1(Az)  for  Az  =  12  km.  Left 
panels  used  multichannel  40  hPa  radiances,  right  panels 
used  channel  71  radiances. 
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and  mountain  wave  effects  are  relatively  more  important 
than  in  the  Antarctic  [e.g.,  Carslaw  et  ah,  1998],  Nonethe¬ 
less,  the  earliest  observations  of  Antarctic  PSCs  and  ozone 
loss  often  showed  both  forming  first  over  the  AP  and 
spreading  out  from  there  \Cariolle  et  al.,  1989;  Ricaud  et 
al.,  1995].  Subsequent  observations  show  PSCs  clustering 
near  the  AP,  yet  Wang  et  al.  [2008]  have  discounted  any 
significant  role  for  AP  mountain  waves  in  this,  stressing 
instead  synoptic  uplift  by  weather  systems.  Synoptic-scale 
dynamics  alone  cannot  explain  the  NAT  outbreak  in  Figure  1 
[Hopfner  et  ah,  2006b],  highlighting  again  the  central  role 
of  AP  mountain  waves  in  initiating  this  particular  NAT 
outbreak. 

[20]  AP  mountain  wave-induced  PSCs  will  also  activate 
chlorine,  leading  to  some  early-winter  ozone  loss  during 
daylight  hours  [e.g.,  Carslaw  et  al,  1998;  Kiihl  et  al, 
2004].  Their  overall  role  in  ozone  loss  is  likely  limited,  of 
course,  since  synoptic  PSC  decks  deep  within  the  vortex 
process  more  chlorine  and  can  account  for  observed  spring¬ 
time  ozone  losses.  However,  NAT  PSCs  are  central  to  the 
denitrification  that  accentuates  this  ozone  loss,  since  STS 
particles  do  not  grow  to  sufficient  sizes  to  settle  out 
gravitationally  to  any  significant  extent  [Tabazadeh  et  al, 
2001].  Thus  the  present  work  suggests  an  important  role  for 
Antarctic  mountain  waves  in  initiating  “NAT  rock”  growth 
and  denitrification  at  freezing-belt  latitudes,  along  similar 
lines  to  that  delineated  from  modeling  studies  in  the  Arctic 
[e.g.,  Mann  et  ah,  2005], 
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